Compton Scattering from the Deuteron and Extracted Neutron Polarizabilities 
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Differential cross sections for Compton scattering from the deuteron were measured at MAX- 
lab for incident photon energies of 55 MeV and 66 MeV at nominal laboratory angles of 45°, 
125°, and 135°. Tagged photons were scattered from liquid deuterium and detected in three Nal 
spectrometers. By comparing the data with theoretical calculations in the framework of a one-boson- 
exchange potential model, the sum and difference of the isospin-averaged nucleon polarizabilities, 
Qjv + Pn = 17.4 ± 3.7 and ajv — Pn ~ 6.4 ± 2.4 (in units of 10~^ fm"^), have been determined. By 
combining the latter with the global-averaged value for a^—fip and using the predictions of the Baldin 
sum rule for the sum of the nucleon polarizabilities, we have obtained values for the neutron electric 
and magnetic polarizabilities of Qf„ = 8.8 ± 2.4(total) ± 3.0(model) and /9„ = 6.5 =F 2.4(total) 
=F 3.0(model), respectively. 

PACS numbers: 25.20.Dc, 13.40.Em, 13.60.Fz, 14.20.Dh 



The electric (a) and magnetic (/?) polarizabilities of 
the nucleon characterize the second-order response of its 
internal structure to applied electric and magnetic fields, 
respectively. Since the polarizabilities manifest them- 
selves in two-photon processes, an excellent method to 
measure them is via Compton scattering experiments. 

The most recent global average (ij for the difference of 
the proton polarizabilities is 

Up- I3p^ 10.5 ± 0.9(stat + syst) ± 0.7(model), (1) 

in units of 10~* fm'^ (which will be used throughout this 
paper) . The sum of the nucleon polarizabilities is usually 
obtained indirectly via the predictions of the Baldin sum 
rule. A recent re-evaluation of this sum rule ^ gives 

ap + /3p = 14.0 ±0.3, (2) 
an + Pn = 15.2 ±0.5. (3) 

Electromagnetic scattering of low-energy neutrons in 
the electric fields of heavy nuclei was used to extract 
an (see Ref. and the references therein). There have 
been many attempts to measure a„ using this method 
0, IE 0i ' but the final conclusion regarding its value 
has remained unclear. Note that this method does not 
constrain /?„ . 

Scattering real photons from neutrons bound in nuclei 
is another way to measure the neutron polarizabilities. 
The use of a deuterium target permits the theoretical 
uncertainties in the interpretation of the experimental 
data to be minimized. Both the reactions — > ^np 
(quasi-free) and ■yd (elastic) can be considered. 

The suggestion to exploit the quasi-free kinematic region 



was made in Refs. An advantage of this approach 

is that it can be carefully tested by comparing measured 
quasi- free proton cross sections with available free-proton 
data. The method has been used in a series of experi- 
ments at MAMI [Hill 113 and SAL 13]. The values for 
the neutron polarizabilities found in these experiments 
are in agreement with each other; however, the sizes of 
the quoted errors differ markedly. The most accurate 
values to date have recently been reported in Ref. |0| as 

an = 12.5± 1.8(stat)tJ;^(syst) ± l.l(model), (4) 
(3n = 2.7t l-8(stat)l?;^(syst) T l-l(model). (5) 

The anti-correlated errors in Eqs. Q and |SJ are due to 
the application of the sum-rule result . In view of the 
model errors contained in Eqs. Q and ij^, confirmation 
of these values is of great importance. 

Elastic photon scattering from the deuteron provides a 
third experimental method for determining the neutron 
polarizabilities. While only the isospin-averaged nucleon 
polarizabilities un — {ap + a„)/2 and Pn = (/3p + (3n)/'^ 
can be measured, this is not a major problem since the 
proton values are rather accurate (see Eqs. ^ and Q). 
Although the first such measurements of elastic -Td scat- 
tering were performed many years ago (see Ref. |l4| and 
the references therein), only two recent experiments at 
Illinois at Ejj^ 49 and 69 MeV and at SAL at E^, 
= 94 MeV [lal have been precise enough to reveal the 
effect of the nucleon polarizabilities. However, the values 
for q;jv — Pn extracted from these measurements are in- 
consistent. With the use of the theoretical model ,2], a 
value of 2.6 ± 1.8 was obtained in Ref. (3| that together 
with Eq. gives q;„ — /?„ ~ — 5.3 ± 3.8. But to describe 
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the data from Ref. [l5j, one needs to increment aM — Pn 
to 7.9 ± 3.8, thus giving = 5.3 ± 7.6. The for- 

mer value is far from theoretical estimates of this quan- 
tity based on dispersion relations, which crudely predict 
that an - I3n ^ ap- (3p (see Refs. [13, EE El)- K also 
contradicts the results obtained from quasi-free neutron 
Compton scattering (see Eqs. (@} and (jSJ). 

In this Letter, we report a new measurement of the dif- 
ferential cross section for deuteron Compton scattering 
performed at MAX-lab. The near-continuous, 95 MeV 
electron beam from the MAX I stretcher ring was used 
to produce tagged photons, in the energy range 50-72 
MeV, with an FWHM energy resolution of ~330 keV 
and a flux of about 3 x 10^ MeV-^s-i H^, EJ. The 
post-bremsstrahlung electrons were momentum analysed 
in a magnetic spectrometer using two 32 scintillator ho- 
doscopes located along the focal plane. They were placed 
such that the central tagged-photon energies were 55 
MeV and 66 MeV. The photon beam was incident upon 
a scattering chamber containing liquid deuterium in a 
cylindrical cell (length 160 mm and diameter 48 mm) [23 | 
made from 125 /im thick Kapton foil. 

Scattered photons were detected in three spectrome- 
ters, each containing a central Nal detector 25.4 cm in di- 
ameter and with depths of either 25.4 or 35.5 cm, placed 
at nominal lab angles of 45°, 125°, and 135° at a distance 
of approximately 0.4 m from the target. The resulting 
solid angle (together with the detection efficiency) was 
precisely determined via GEANT simulations. The en- 
ergy resolution of the Nal spectrometers ranged from 6 
to 8% at a photon energy of 60 MeV. The gain stabil- 
ity of the Nal detectors was continuously monitored and 
instabilities were corrected for using a Light Emitting 
Diode (LED) system j2^. In turn, the stability of the 
LED system was monitored and verified by examining 
the location of the cosmic peak on a run- by-run basis. 
The data acquisition was started by an event in any one 
of the Nal detectors, which provided gates for charge in- 
tegrating ADCs (Analog to Digital Converter) and start 
signals for TDCs (Time to Digital Converter) used for 
time-of-fiight (TOF) measurements. The stop signals for 
the TDCs came from the focal-plane detectors. In order 
to monitor the number of pile-up events in the Nal crys- 
tals, a 250 MHz Flash ADC was used. Less than 1% of 
the events were affected by pile-up. 

Data were collected over 8 weeks, divided typically into 
two-week run periods. For each run period, tagging ef- 
ficiencies (~20%) were measured using a Pb/SCIFI de- 
tector [2^. In addition, the responses of the Nal spec- 
trometers were measured by placing them directly in the 
reduced intensity photon beam. For each detector, the 
aforementioned GEANT simulations were phenomeno- 
logically broadened to match the measured responses. 
This broadening was then folded into second-stage simu- 
lations of the in situ responses and solid angles 0| . A 
typical result is shown by the dashed line in Fig. ^ 
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FIG. 1: Smiss spectrum at O^^'^ = 126° and E-^ = 55 MeV. 
The two vertical lines indicate the region of interest (ROI) 
used to determine the yield. The solid line represents the 
sum of the fitted responses (see text). 



This missing energy (-Bmiss) spectrum was obtained by 
summing over a 10 MeV interval centered at 55 MeV. 
-Bniiss was defined as the difference between the tagged- 
photon energy (corrected for the Compton scattering en- 
ergy shift) and the energy registered in the Nal detec- 
tor. The subtraction of random events was performed 
using two independent methods: one used the TOF be- 
tween the scattered photons and the post-bremsstrahlung 
electrons, while the other employed the non-physical 
(i?iniss < 0) region of the energy spectrum. An aver- 
age normalization factor was employed which resulted 
in a 4% systematic uncertainty in the extracted cross 
sections. The experimental data were fitted within the 
region of interest (ROI) considering the contributions of 
elastically scattered photons from liquid deuterium, from 
the Kapton cell and the scattering chamber windows, and 
from ice (H2O) which built up on the cell during the run 
period. The contribution from the cell and windows was 
calculated using known differential cross sections for car- 
bon and oxygen [25I I2H . W\ . The contribution from ice 
was fitted using data from the largest detector at the 
backward angles and extrapolated to the other detector 
angles. The careful quantification of the experimental 
background due to the empty target and ice resulted in 
an average improvement in the of the fit functions 
of 20% for the backward detector angles. The observed 
contributions correspond to ice-layer thicknesses of about 
100 /xm per cell endcap. The contribution from inelastic 
scattering (which begins at -Bmiss = 2.2 MeV) was kept to 
less than 3% via the narrow ROI. The systematic uncer- 
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TABLE I: CM differential cross sections for deuteron Comp- 
ton scattering. The first uncertainty is statistical and the 
second is systematic. 
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tainties for the experiment (~8%) arise from the tagging 
efficiency (5%), the product of the soUd angle and the 
detection efficiency (4%), the random background sub- 
traction (4%), the contamination of the ROI by inelastic 
photon scattering (3%), and the target thickness (2%). 

Results for the center-of-mass (CM) cross sections are 
given in Table |l| and displayed in Fig. [3 Earlier results 
from Illinois [l5j and SAL [l6j are also shown. Since the 
Illinois data were obtained at slightly different photon 
energies (49 MeV and 69 MeV), we extrapolated them 
to our energies using the theoretical model presented in 
Ref. "^I. Good agreement with the Illinois data is clearly 
demonstrated. 

The polarizabilities of the nucleon have been deter- 
mined using the theoretical model presented in Ref. 0- 
In this model, apart from effects due to the one-body 
interaction (which contain the effects due to the polar- 
izabilities), two-body effects due to the NN interaction 
and related Meson-Exchange Currents (MEC) have been 
included via a series of one-boson exchanges which consti- 
tute the Bonn potential. After making a two-parameter 
fit to the data in Tabled with the use of this model, we 
obtain for the sum and difference of the isospin-averaged 
nucleon polarizabilities 



17.4±3.7(stat 
6.4 ± 2.4(stat 



^syst), 
-syst). 



(6) 
(7) 



with x^/^doi — 7.5/(18 — 2). The quoted uncertainties 
are the statistical and systematic uncertainties taken in 
quadrature. The obtained sum un + Pn is in agreement 
with the values given in Eqs. (0) and This indicates 
that the systematic uncertainties (including the model 
uncertainty) are well understood. 

Having extracted un and /^at, we invoked the accurate 
proton values for Up and j3p to derive the neutron polariz- 
abilities. In doing this, we relied on the precise sum-rule 
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FIG. 2: CM differential cross sections for Compton scatter- 
ing from the deuteron. Filled circles: pres ent experiment; 
open circles: extrapolated Illinois results llSll (see text); open 
squares: SAL results [l^ . The error bars represent the 
quadratic sum of the statistical and systematic uncertain- 
ties. Solid and dashed curves are the calculations of Ref. 0| 
(dashed: ajv — Pn ~ 0; solid an = 11.9, /3jv — 5.5, see 
Eqs. (|SJ, Q). Dotted and dashed-dotted curves are the pre- 
dictions of Refs. 28] and ^29j . respectively, extrapolated to 
Qfiv = 11.9, (3n = 5.5. 



prediction rather than on our result (jnj. Making use 
of Eqs. CJ -® and 0, we obtain 



an = 8.8 ± 2.4(stat + syst), 
/?„ = 6.5 T2.4(stat-K syst). 



(8) 
(9) 



If all currently available data (Illinois, SAL, and Lund) 
are fit using the theoretical model of Ref. 2J , the follow- 
ing "global values" may be inferred: 



aN + Pn = 16.7± 1.6(tot), 
ctN - Pn 4.8± 2.0(tot), 



(10) 
(11) 



with xV^dof = 38/(29 - 
rule and Eq. (P), 

a-n = 

Pn = 



2). Thus, using the Baldin sum 



7.2±2.1(tot), 
8.lT2.1(tot). 



(12) 
(13) 
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Here both statistical and systematic uncertainties have 
been combined, the latter being taken into account 
through a rescaling of measured cross sections within 
their normalization uncertainties. 

Model uncertainties in the extracted values of the po- 
larizabilities can be partly understood by comparing dif- 
ferent calculations of da/df^ 0,|ii|2i|3a|3ll'3^. Fig.H 
shows the three most recent predictions |2, 12^ l2Bj at 
fixed aN — 11-9 and = 5.5 as given by (0) and 
(0. At the energies and angles of the present experi- 
ment, there is reasonable agreement between the poten- 
tial model of Ref. |3] and the Of^f) Chiral Perturbation 
Theory (ChPT) results of Ref. j^l. At the SAL energy 
and angles, this agreement is poorer. Agreement with 
the 0{p^) ChPT calculation of Ref. jH is worse. 

Within the framework of the potential model used 
here, the main uncertainties which contribute to the ex- 
traction of the polarizabilities arise from evaluations of 
the MEC and seagull terms. The majority of these uncer- 
tainties may be collected into the Thomas-Reiche-Kuhn 
sum rule enhancement parameter k which determines the 
magnitude of the seagull contribution 2]. Depending on 
the NN potential used (Paris, Bonn, OBEPR, Argonne 
vl8, Nijm93, or CD-Bonn) the value of k varies from 
0.44 to 0.51. This leads to variations in the differential 
cross section of about 4% which result in variations in 
the extracted value of ajv of about ±1.5. We adopt this 
as an estimate of the model uncertainty in a^v . Accord- 
ingly, the model uncertainty in the derived value of 
is about ±3. The model uncertainty in /?„ is anticipated 
to be smaller than that in a„, as the changes in k given 
above affect (3m to a smaller degree. 

The values obtained for q;„ and /3„ are in reasonable 
agreement with those found in quasi-free Compton scat- 
tering from the neutron (see Eqs. and ©). How- 
ever, we observed the tendency of the technique of elas- 
tic Compton scattering from the dcutcron to give smaller 
values for — /3„ than did the technique of quasi-free 
scatterin g. P redictions for q:„ in the framework of ChPT 
at 0{p^) 33] , the heavy baryon ChPT model with the A- 
isobar included [SJLand the so-called "covariant dressed 
K-matrix model" [35j give values of 13.0 ± 1.5, 16.4, 
and 12.7, respectively. All are somewhat larger than our 
value. The situation concerning /3„ is less clear. Reson- 
able agreement with ChPT results of 7.8 ± 3.6 and 
9.1 |3J| is demonstrated. Ref. [s^ suggests a noticably 
smaller value of 1.8. 

In summary, differential cross sections for deuteron 
Compton scattering have been measured at MAX-lab. 
The data were used to extract neutron polarizabilities in 
a model-dependent analysis. The extracted values for Q!„ 
and /3„ are consistent with those obtained from the Illi- 
nois data [lB |: however, they are inconsistent with those 
resulting from an analysis of the higher energy SAL data 
|l6|. The present analysis thus confirms the previous 
observations 0, |i32j that the available models cannot rec- 



oncile data obtained at photon energies of about 60 MeV 
with those obtained at about 100 MeV. 
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